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SUMMARY 


A method has heen developed to determine the flow field of a 
body of revolution in steady separated flow (e.g. at high angle of 
attack) , including the prediction of the normal force and pitching 
moment . 

The method makes use of a computer to integrate the various 
solutions and solution properties of the sub-flow fields which make 
up the entire flow field. Thus, a finite difference solution to the 
complete Navier-Stokes equations is not employed. 

The method utilizes two ideas in this approach to a steady three- 
dimensional separated flow: the unsteady cross flow analogy which 

reduces the given three-dimensional steady separated flow to a two- 
dimensional unsteady separated flow and a new solution technique for 
the latter problem. This technique employs a wake description of 
discrete point vortices arising from the separation of shear layers at 
the surface. The point vortices convect aua diffuse downstream to form 
an unsteady wake. Thus the mathematical model follows directly from 
the physics of the wake evolution with time. 

The overall technique is applied, employing a computer, to two 
test cases: an ogive-cylinder and an ellipsoid of revolution at low 

speeds. Force and moment data are obtained and are found to agree well 
with experimental data, particularly at high angles of attack where 
inviscid theory is invalid. Separation regions and wake patterns are 
found which agree with available experimental findings. 


1 . 0 Introduction 

The objective of the present work is to develop a numerical 
computation method to determine the flow field of a body of revolution 
in steady separated flow (e.g., a body at high angle of attack), 
including the prediction of forces and moments. 

There are many computational techniques of varying accuracy for 
the flow about a body with no separation such that the assumptions of 
inviscid flow and linearized equations of motion can be utilized 
(references 1 & 2, for example). For such flow fields the techniques 
are based upon the availability of fundamental solutions of the 
equations of motion in integral form but with unknown singularity 
distributions. The computer is then used to determine specific 
solutions for the particular bodies which meet imposed boundary 
conditions and to evaluate the resulting numerical solution to yield 
pressure distributions, forces and moments. For flow fields in 
which separation must be taken into account, the assumptions of 
inviscid flow and linearized theory are invalid and the usual approach 
is to solve the Navier-Stokes equations in finite-difference form by 
means of a computer (references 3, A, 5). At present, the solutions 
for inviscid, linear flow are more developed and hence, more reliable, 
and require far less computation than the full finite difference 
solutions. 

It is intended here to adopt the former approach; i.e., 
computer implementation of known solutions and solution properties. 

This approach is believed to promote an increased physical understanding 
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of the problem, to reduce the role of the computer as a black box, 
and to minimize the computational time and storage. 

The obstacle in pursuing such a course is, obviously, the 
relatively very few solutions and solution properties for viscous 
flows compared to the situation for inviscid flows. To overcome 
this obstacle, two techniques are employed. 

The first is the unsteady cross flow analogy (or viscous cross 
flow analogy or impulse flow analogy) , which is based upon the 
assumption of an equivalence between three dimensional steady flow 
and two dimensional unsteady flow. As explained in more detail 
below, this assumption is a heuristic one based primarily upon 
experimental data. With this technique, the three dimensional 
steady separated flow problem is reduced to a two dimensional unsteady 
separated flow problem. 

The second technique, applicable to the latter problem, is 
based upon the assumption that the two dimensional unsteady wake can 
be described by a distribution of inviscid point vortices superimposed 
on the unseparated potential flow solution, suitably modified by 
diffusive effects. The argument for this is based upon the time 
evolution of the wake. 

Part I of the report treats the basic theory, applications and 
results while Part II describes the computer program. 

In the following sections of this part, Part I, the basic theory 
consists of the description of the unsteady cross flow analogy and the 
solution technique for the flow in the cross plane, a two dimensional 


unsteady flow field. The theory is applied to two geometries, an 
ogive-cylinder and an ellipsoid of revolution and results for normal 
force, pitching moment, vortex-wake patterns and separation regions 
are compared with experimental data. 
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2,0 Unsteady Cross Flow 


2.1 General Concepts 

The problem is to determine the flow field and to predict the 
forces and moments on a three-dimensional body of angle of attack a 
in a steady uniform flow (Figure 1) at high Reynolds number. The 
angle of attack a is of arbitrary magnitude such that the flow may 
be separated. Thus linearized theory is not sufficient. 

„For aerodynamic bodies at small angles of attack, the three 
dimensional boundary layer adheres to the entire body except in the 
neighborhood of a small base. For this situation the boundary layer 
may be ignored allowing for an inviscid flow approximation and, for 
slender bodies, application of linearized theory. For higher angles 
of attack (references 6-10) the boundary layer separates and may 
separate to such an extent that flow separation exists over a substantial 
portion of the body surface. In this situation, inviscid theory is 
no longer applicable (even though the body is slender) and there is 
no viscous theory to describe this complicated three dimensional flow 
field. 

This problem is approached herein by using the unsteady cross 
flow analogy applied to a body of revolution (references 11-23) . The 
essence of this technique is found in applying the (dimensional) 
transformation [x*, y*, z*, £*] ■ [x*, y*, z* + Wt*, t*] and considering 
the flow field in the plane z* ■ 0 (see Figure 1). In this plane, 
the flow field is that induced by a circular cylinder of time-varying 
radius (correlated with the thickness distribution of the body) in a 



uniform, viscous, incompressible free stream of velocity U ■ V sin a 


Information available from the relatively simpler two-dimensional 
unsteady flow is thus applicable to the three-dimensional steady flow 
field. In particular, the two-dimensional unsteady drag distribution 
with time is applicable. This can be transformed to the longitudinal 
normal force distribution of the three dimensional steady flow field 
yielding the normal force and pitching moment on the given body at 
angle of attack. 

The correlation of the two dimensional unsteady results and the 
three dimensional steady results lacks a theoretical foundation 
(references 13, p. 415 and 24, p. 464). Basically, there exists no 
general transformation between the two dimensional unsteady and the 
three dimensional steady Navier-Stokes equations. There exist parti- 
cular transformations such as in slender body theory and its extension 
to hypersonic flow, the equivalence principle, and it is of note that 
these are inviscid theories (shocks entering by jump conditions). So, 
to this date, the prime argument for the validity of the unsteady cross 
flow analogy lies in the similarity of the flow fields obtained, 
usually by vortex visualization techniques, and the valid results 
achieved, although the latter can be argued against on the grounds of 
empiricism. 

The question of the existence of a transformation between three 
dimensional steady and two dimensional unsteady Navier-Stokes equations 
can be considered in terms of a basic property of the unsteady cross 
flow analogy, namely, the absence of upstream effects; the force 
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A 

distribution on the three dimensional body over the length (0, z*) 
is independent of the force distribution over the length (z*, Si ) . 

For this to be true the steady Navier-Stokes equations would have 
to be parabolic or hyperbolic with z* as a time-like coordinate. 

If the three dimensional steady flow field is decomposed into 
the outer inviscid flow, the boundary layer, and the wake, the 
governing equation would be the Euler equations, the boundary layer 
equations, and the complete Navier-Stokes equations, respectively. 

For incompressible flow (reference 25) , the outer flow is governed 
by elliptic equations, the boundary layer by parabolic equations with 
z time-like, and the wake by elliptic equations. Thus only the boundary 
layer region is compatible with the unsteady cross flow property. 

For supersonic flow the theory is not as well developed. The 
outer flow is governed by elliptic equations (which mathematically 
may be ignored for a sharp pointed body) and hyperbolic equations. The 
boundary layer is still parabolic but the complete Navier-Stokes 
equations in the wake still have an elliptic property. 

Thus, aside from shocks, the requirements for the validity of 
the unsteady cross flow analogy are partially met in supersonic flow 
in that the outer flow is for the most part hyperbolic but in the 
dominant region, the wake, the requirements are not satisfied. 

While the boundary layer of the three dimensional steady flow 
field is parabolic, there is the conflict between three dimensional 
and two dimensional boundary layer states (reference 13, p. 415). For 
usual three dimensional flow fields, the three dimensional boundary 
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layer is for the most part turbulent. Yet the pertinent Reynolds' 
number for two-dimensional unsteady flow, a function of the Reynolds' 
number of the three dimensional steady flow, could be such that a 
two dimensional laminar boundary layer theoretically exists. This 
mathematical result has no physical interpretation. For the particular 
three dimensional steady flow fields considered herein and for the 
particular equivalence between three dimensional steady and two-dimensional 
unsteady flow fields used herein, the pertinent Reynolds' number 
determined for the latter flow field are such that laminar boundary 
layers theoretically exist. Thus the use of laminar two dimensional 
unsteady boundary layers are employed based upon this formal justification. 

In addition there is the discrepancy between the concept of 
two dimensional and three dimensional separation of boundary layers, 
further complicated by the unsteadiness of the two dimensional boundary 
layer . 

There exists a body of work employing the application of the 
unsteady cross flow analogy to a delta wing of low aspect ratio. This 
problem has the advantage in that the sub-problem of the two dimensional 
unsteady flow (about a flat plate or thin ellipse as opposed to a 
circular cylinder) is simpler since the location of the separation 
point is known a priori (i.e. the flow separates almost immediately 
at the edge). In addition, the delta wing can be approached by 
assuming a conical flow in the three dimensional steady case which 
bears resemblance to the two dimensional unsteady case (reference 26) . 

This area has received a great deal of attention (references 27-30) 
but the present work cannot make use of these simplifications. 



The work of Hill for bodies of revolution (reference 31) , while 
not employing the explicit cross flow analogy, bears upon the problem 
of separated flow. The theory is completely inviscid and so the 
approach is a limited one. 

In this context then it is not surprising that the correlation 
between the two flow fields has been achieved in various ways, often 
requiring ad hoc experimental data. The correlation used in this 
work (see Section 2.2) is the direct one of equating two dimensional 
unsteady drag (pressure and viscous, although the latter is usually 
negligible) to the normal force distribution. 

The sub-problem of two-dimensional unsteady flow while simpler to 
treat than the original three dimensional steady one, is still a 
difficult problem since a boundary layer forms, separates, and a wake 
is formed about a cylinder of varying radius. This sub-problem, discussed 
more fully in Section 3.0, is treated with a new approach compatible 
with the originally defined approach of using the computer with known 
solutions. With the possible exception of reference 17, (which uses 
a finite difference approximation to the Navier-Stokes equation) this 
approach yields the most complete model (for high Reynolds' number) in 
that viscosity enters directly and empiricism is minimized. 

2.2 Transformations 

In the absence of a transformation between the three dimensional 
steady and two dimensional unsteady Navier-Stokes equations, the 
unsteady cross flow analogy is applied by utilizing simple transfor- 
mations of the geometries and the forces. 



The geometry transformations obtained as follows: given the 


three dimensional steady flow field of Figure 1. If a velocity 
W = V cos a is superimposed on the steady three dimensional flow 
field of Figure 1, the flow field of Figure 2a results. Then if 
attention is confined to the plane z* = 0, the situation is that of 
an unsteady two dimensional flow field with a constant velocity at 
infinity and body radius changing with time, Figure 2b. 

The geometry transformations are then, with ( )* as a dimensional 
quantity, 

A 

a*(t*) - r o *(z*), t* = ~ J 0 < < A, 0 < t* < £ (1) 

for arbitrary r Q *(z*). 

Now r Q * is normalized by d/2. In the unsteady two dimensional 

problem, a*(t*) must be normalized- by a length characteristic of the 

problem. Let this characteristic length be a to be defined for each 

fWM 

particular body but on the order of (W/i) a*(t*)dt* ■ (1/i.) • 

tl JQ 

r *(£*)d£*. 

Jo ° 

Letting t = Ut*/a, the above equations in nondimensional form 

are 


fi (t) - (f^) r o (z), t 


/ d f \ a ^ ^ a ^ ^ ^ ^ / df fan u 

(-=— tan a) z; 0 - z _ 1, 0 : t : » 

a a 


The force transformation is obtained as follows: given the 

three dimensional steady flow field of Figure 1. The normal, force, 
N, and the pitching moment, H, are to be determined utilizing the 
unsteady cross flow analogy. 
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Consider the three dimensional steady flow field. Let N(z*) 
be the force normal to axis of symmetry on the section of the body 


from the nose to z* such that N(£) = N. Then 


N = 


™di* 

dz* 


( 3 ) 


The moment about the nose is then, positive nose up, 

rSL dr 


M = - 


[z* + r * -prj-] prr dz* 
o dz* dz* 


(4) 


and the moment about an arbitrary point [x,y,z] = [0,0,X&], 0 £ X - 1» 
is 


M - M + X£N 

A 


To reduce these to coefficient form let 

„ dN 


(5) 


= n <2> 


N 


"M.X 


dz* 

0 

( 6 ) 

pV 2 d 


2N 
- 2 

(7) 

pV S 


2M. 


A 

( 8 ) 


where d is the maximum diameter, l is the length and S ■ ird /4 the 
frontal area. 

Normalizing z* by A, equations (3), (4) and (5) become 


, 4f f 1 

"N tt J C n 

J n 


dz 


Si, x 


4f 

TT 


1 r dr 

[z + — 5 - c dz + XC 

o 4f 2 dz n N 


(9) 


( 10 ) 
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The above coefficients are found by solving the two dimensional 
unsteady problem. Within this problem, the pertinent force is the 
drag as a function of time. (It is assumed that the times are not 
long enough for any oscillating lift to develop). 

The drag is given by 


D(t*) 


2tt 


p*cos 6 a*(t*)d0 + 


2tt 


t*sin 9 a*(t*)d6 (11) 


With a coefficient of pressure 

2(p* - P* ) 


PU 


the coefficient of drag 


C D = 


2D 


pU 2 (2a) 


is obtained as 


( 12 ) 


(13) 


C D (t) = 2 


2tt 


C cos 0 a(t)d0 + — 
P 1 


2tt 


t sin 0 a(t)d0 


(14) 


2 

where x * 2x*/pU . 

| 

Then the basic assumption in the force transformation is 


j£*(z*) = D(t*), z* = Wt*; v = v,p = p, p oo = p oo (15) 

l 

i 

or in non-dimensional form 


2a .2 , v 

c ( 2 ) = — sin a C (t) ; z 
n a u 


at 


df tan a 


(16) 


with Re 2DUS = Re 3Ds sina and the same free stream static values. 
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This assumption ignores longitudinal shear effects. But for 

aerodynamics shapes the fineness ratio f ■ £/d is usually large, thus 

reducing the contribution of longitudinal effects. 

To further check the assumptions, comparisons were made between 

the experimental pressure distributions of the three dimensional 

steady flow and the theoretically derived two dimensional unsteady 

ones. To make this comparison let N (z*) be the normal force due 

P 

to pressure only on the segment of the body from the nose to z* in the 
three dimensional steady flow. Then 


dN r 2ir 

E „ 

dz* 


p*cos 6 r Q *(z*) d9 


(17) 


In terms of the coefficient of pressure 

2(P* - P J 


Cp 


and the sectional coefficient 


r \ x .2 

pV 


dN 

2 — E - 


(18) 


dz* 

n,P j$V 2 d 


(19) 


c (z) - j 
n,p 2 


2ir 


C (9,z) cos 0 r (z) d0 
p o 


( 20 ) 


In two dimensional unsteady flow, let the portion of D(t) due 
to pressure be 


D p (t) 


'2tt 

p*cos 0 a*(t*) d0 
' o 


( 21 ) 


In terms of the coefficient of pressure as in equation (12) and the 


coefficient of pressure drag 


D,p 


2D(t*) 

pU 2 (2a) 


"D,p 


1 

2 


2ir 


C cos 9 a(t) d0 


( 22 ) 


(23) 


From equation (2) and the pressure contribution to equation (16) , 


C (6,z) = C (9,t) sin^a 

P P 


(24) 


However, as is seen in the following section, this comparison 
can only be made qualitatively. 


3.0 Two Dimensional Unsteady Separated Flow 
3 . 1 General Concepts 

The problem is to find the flow field induced by a circular 
cylinder of time-varying radius in a uniform stream of a viscous 
incompressible fluid (Figure 2b) . The results are to be applied to 
predict the forces and moments of a body in a three dimensional 
steady flow field. 

The flow field is that of a uniform flow with no body present 
for t<0. At t ■ 0 a circular body appears at the origin with a 
radius a(t) such that a(0) - 0. During the initial stages a boundary 
layer (assumed to be laminar) is formed on the circular cylinder. 

At later times the boundary layer separates from the cylinder, 
creating vorticity in the flow field immediately after the cylinder. 
Thus a wake is formed. 


17 



In previous work employing the unsteady cross flow with a body 
of revolution (references 11-23) , the two dimensional unsteady 
solution was derived from a direct use of experimental data or the 
use of inviscid solutions with ad hoc experimental data to describe 
the viscous effects. Of note is the work of reference 17 which used 
the finite difference approach for the two dimensional unsteady 
problem. 

The approach in the present work is unique in that the viscous 
element, the primary mechanism being the separation of an unsteady 
boundary layer, is included without resorting to a finite difference 
approximation to the complete Navier-Stokes equations nor to ad hoc 
experimental data. This approach can be described as follows: 

The problem at early times is sometimes referred to as the cone 
problem for which no solution exists. It will be approximated by the 
small time analytical solution (since numerical integration follows, 
an analytical solution is required for initial conditions) for a 
circular cylinder of constant radius impulsively started from rest in 
a viscous incompressible fluid (reference 32) . Since the latter can 
be achieved experimentally by plunging a circular rod into a flowing 
stream, this approximation is valid for a blunt nosed three dimensional 
body. 

At later times, the solution is obtained by using the technique 
employed in other work in obtaining the solution for large times for 
a circular cylinder of constant radius impulsively started from rest 
(this work will be referred to as "the constant-radius-solution") , 
and modifying it for a time-varying radius. 



The problem for a constant radius cylinder which has received a 

i 

great amount of attention, is considered a classical problem in 
separated flow (reference 33). There have been two primary approaches: 

| a solution to the Navier-Stokes equation in finite difference form 
(reference 34) and a purely inviscid solution (references 14,15). The 
i latter while simpler still requires a computer and must employ 

I 

| empiricism for the separation phenomenon, a viscous action. 

The approach adopted in the constant-radius-solution was to 
use the computer with known solutions (and matching different solutions 
based upon the physics of the flow field) while retaining viscous 
phenomena. This approach has succeeded except for the case of the 
two dimensional unsteady boundary layer where a finite-difference 

solution was required. 

I 

(An approach midway between a finite difference approximation 
to the Navier-Stokes equation and the model used here is the work of 
Wu (reference 5) (which did not use unsteady cross flow concepts). The 
difference is that Wu solved the Navier-Stokes equations in the wake 
alone while here a set of point vortices was used. Wu's model, while 
better theoretically, has not been successful at the values of Reynold* 
I numbers used here). 

1 The particulars of the approach of the constant-radius-solution 

are as follows (see Figure 4) : For the time period during which the 

cylinder moves about a third of its radius the boundary layer is 
fully attached with classical potential flow outside. Starting at 
some time before separation, the flow field is given as the 


unseparated time dependent boundary layer solution (reference 32) , 
evaluated at this time, as driven by the classical potential flow 
about the cylinder. The solutions for the laminar boundary layer 
for later times are then obtained by numerical integration of the 
boundary layer equations employing Hall's scheme (reference 35). 

With succeeding times, the boundary layer separates starting at the 
rear stagnation point and the point of separation moves upstream. 

At the point of separation, the vorticity flux across the boundary 
layer is found and during a time At, corresponding to the time of 
integration of the boundary layer, this flux is summed up and 
represented by a point vortex. This point vortex is convected 
downstream and its effect is modified by diffusion. Thus, the wake 
is formed of these point vortices superimposed on the original 
potential flow. The wake grows, being continually fed with vorticity 
by the separating boundary layer, while outside of the wake and the 
boundary layer, irrotational flow exists. This situation continues 
until the steady state region is attained. Physically a rigorous 
steady state is never achieved due to instabilities in the wake and 
mathematically this may only be attained in infinite time. The work 
of the constant-radius-solution reveals the neighborhood of the steady 
state is achieved in a finite time (Figure 5) but for its present 
use, 0 < t < (df/a) tan a, the steady state region does not enter. 

It is to be noted that the no-slip condition is satisfied on 
the portion of the cylinder covered by the boundary layer but not on 
the remaining portion. The wake model as depicted above while 


reducing the surface speed of potential flow is not sufficient to 
yield a surface speed sufficiently less than free stream speed to 
approximate to a no-slip condition. Thus there is postulated a shear 
layer on the rear of the cylinder, similar to a boundary layer but 
with a rotational outer flow, to satisfy no-slip. This rear shear 
layer may separate, creating vorticity in the wake of the sign 
opposite to that shed by the boundary layer separation. 

For the present problem, the above physical situation is modified 
by the presence of a time varying radius a(t). 

With this approach the flow field is decomposed into a set of 
sub-flow-fields and their interactions. The governing equations and 
their solution are given in the following. 


3.2 Outer Flow 

The flow outside the boundary layer rear shear layer (see Figure 4) 
is assumed to be the classical potential flow (uniform flow and source 
and doublet at origin) over a cylinder plus the potential flow induced 
by a set of point vortices. 

In terms of the non-dimensional quantities [x,y,t] = [x*/a,y*/a,Ut*/a] , 
r = r */a*, $ + i 'f = ($* + iY*) / Ua, to - a co*/U, r = r*/Ua, the 
equations for the outer flow are 


DE 

A'?(x,y,t) = u)(x,y,t) 


IC 

t 

II 

o 

-6 

II 

(25) 

BC 

1) 

r = a : T = -aa6 



2) 

r oo ; m y 
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where w(x,y,t) = £ f^ k 6(x-x^ k ) 5 (y _ yg k ) 


(6 being the delta function) such that 


wdS " l r Bk* 


The w(x,y,t) arises from vorticity created by the separation of the 
boundary and rear shear layers. 

The solution is 


w = $ + iy 


2 

z + - aa log z + i l 


Bk 

2tt 


log 


< z iVi 


2 

a z 


z - 


Bk 


'Bk 


'Bk 1 


where 


z = x + iy = re 


ie 


10 


V = x 8k + iy 8k ‘ V 


6k 


10 


z - z 6 k ’ + 1(y " y Bk ) " r lSk e 


16k 


2 

a z 


z - 


Bk 


2 

a x 


= x - 


Bk 


+ i 


a y 


y - 


Bk 


i9 


r 2Bk® 


2$k 


'Bk 1 


Bk 


Bk 


such that when 


l z l ‘ r 2Bk ’ a WV and * * °- 


Now 

dw ~ O i 1~0 “i0 / o i j o\ 

— = -u + iv = e (-v + iu ) 

dz 

i 

The form of the equations used in the computations, as in Part II, 
obtained as follows; since 


(26) 


(27) 

are 
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4s - i -i T -S&+ i l 

dz 0 9 L 


F gk f 1 

2 ’ l z-z 6k 
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then 


- 0* + «• - 1 - a 2 t^-y 2 - 21 ^ ] 
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(x-iy) 
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2 
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2 
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'2gk 
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(28) 


and 


o , • o a 

-v + iu = e 


2 -16 • .. r_. ( „~ i6 iBk ‘ ie 2Sk -18 'l 

-V- - — + le I ( 5 + — 

r 2 1 6 2 * l r 16k r 2Bk r J 


(29) 


Since 
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-v° + iu° - e 


ie 


2 -ie • r . f r 

a e aa . j v 8k J _ 
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- — + i l 
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2 

r 18k 


r - f- e^'V 


'8k 


2gk 
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such that when r * a. 


-v£ + iuS * -a + i 2 sin 6 + i £ 


„ . 2 2 

r 8k J a "V 


„ 2n \ 2 

3 ' a r 18k 


+ - 


Non-dimens ionalizing the pressure as p ■ 2p*/pU , the Bernoulli 
equation is 

(uS + vS ) 


- *, t + 


+ f - f(t) 


( 


Then since 


* ■ cos 6(r + - “ 1° 8 r - jj iir (8 isk - e 2 *k + 8 - 9 sk 

P 


*’t - cos 6 - |j- (aa) log r - I jSL 1$^ - - e^} 
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where 
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18k 


18 k 


f d x " a2x 


ftv a r a2x «v'i y- a2 y ft v^ . i 



Since 




gk dt 


(aa) lim log r. 


p P 

f(t) ■ ■ I 2? V + 2 + T + dF <aa) lin 108 r 


such that the pressure coefficient on the surface 


C h p - Poo = [2$, - (uS 2 + vS 2 ) + 1] 

P ^ r*=n 


- 2 £ *r V + 2 If (aa) llm 108 r 


Considering only those terms contributing to drag 


C p - 4 i cos e - 2 l ^ [e 16k - e 2sk ) r=a 

3 


- 2 sin 0 


2 2 

T at a -it. , 

+ y _§k gk l 

4 2tt 2 2 / , + a 

6 _ a r iek () 


Note that the term in equation (34), lim log r, precludes an absolute 
determination of C _ 


3 » 3 Boundary Layer 

The unsteady two dimensional Navier-Stokes equations in polar 
coordinates are (nondimensionalizing as in the outer flow and with 


x = 2x*/pU ) , 




DE 


(rv), r 
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where 


with 


Du . uv p, $ . 2 , . , . 2v, 0. 

Dt + ^ 2^ + 5 [L[ul + — 1 


Dv u 2 p ’r . ,2., , , 2u, e. 

Dt ~ r — + S Il[vl 1 


n u ( )>a 

Dt " < )j t + r ” + v< >•, 
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26 2 u, 


and where 


1 

Re 2DUS 



(36a) 


IC t ■ t > 0; u = f(r,9), v = g(r,0) (36b) 

o 

and boundary conditions at the surface 

BC r = a(t): u = 0, v = a (36c) 

Defining the inner variables 

(v - — ) 

[u,v] = [u, } , [e,T,l] = [0, , t] (37) 

substituting into equations (36) and letting {+0, one obtains the 
(laminar) boundary layer equations. 
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DE 


u ’e 


+ v,- = o 
a ’r 


- - - uu - 9 aru,- • _ p,- 

“•H + vu * r + -j r i + t u --2r + ;: ’ 


rr 


with 




" P, r ’ “ " ~r • T = 2 6 


IC 


t = t : u = f ' (r,9) , V = g’ (r,0) 


BC 1) r = 0: u=v=0 

2) r -> «: u ■> u£(e,t) 

BC2) being derived from the matching concept. 

As such the Hall scheme of integration (reference 38) is not 
applicable. If the transformation is made, 

i 

u = au 


i - a - 
V = V - — r 
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then the problem is 


DE 


u 


’0 , i a 

— + V = 
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IC t Q - t u : u 1 = f"(r,e), v 1 = g"(r,0) 


(38a) 


(38b) 


(38c) 


(39) 


(40a) 


(40b) 
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BC 1) r - 0: u 1 - v 1 - 0 

2) r -*• u 1 -+ auJJ(9,t) 

where p,- has been evaluated at r -► » by virtue of p,- « 0 (40c) 

a r 

Except for the source term "-a/a" in the continuity equation and 
2 

the factor "1/a " multiplying the ( ),- terms, the form of the equation 

is identical with the classical equations. Thus modifying the Hall 

method by the inclusion of a source term and modifying finite 

2 

difference in 6 terms by the factor 1/a (t), the Hall scheme is applicable. 

3.4 Boundary Layer Separation Region 

The numerical integration scheme used to obtain solutions to the 
unsteady two dimensional laminar boundary layer equations, being based 
upon boundary layer approximations, is not valid at the separation 
region. In addition, a definition of the separation region in unsteady 
flow is a controversial issue (references 36-38) . The approach used 
here is the classical one of carrying over the definition of separation 
for steady flow, i.e. t = 0, to unsteady flow. 

Employing x ■ 0 as the definition of separation, the numerical 
scheme used did predict a meaningful separation region (based upon the 
constant-radius solution) . In the scheme actually used the particular 
grid point on the cylinder immediately upstream of the point at 
which x f 0 was used to define the separation point, 6 g . 

3.5 Boundary Layer - Outer Flow Interface 

Let the separation point of the boundary layer be O^t^) = 9^ 
where t^ are the discrete times of the numerical integration. Physically, 
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at this point, the vorticity created in the boundary layer passes into 

the outer flow. This is found as follows: 

The vorticity flux across the section of the boundary layer at 

0 , is 
sk 


r -t<V - 


w u dr 


( 41 ) 


Since <u = u,-/6, u = u and dr = <5dr - in the boundary layer 

.»2 


r ’t = 


-2 

_ _ u 
u,- udr= F 


u 0 


( 42 ) 


using no-slip. 

During the period t^ < t < t_^ , this vortex flux is summed 
into one point vortex of strength. 


r °i' l_ r, t * fc k* At k 5 At k = C k+1 " \ 
6=k 


( 43 ) 


where S denotes the point vortex label. 

Considering this vortex to be in the outer potential flow, its 
images are born at the same time. The surface velocity due to this 
just born point vortex and its images is ~ r gj c / irm j c where m^ is its 
distance from the surface of cylinder into the flow. If this velocity 
is used to cancel the velocity at the surface induced by the outer 
flow minus this new born vortex, then its position is 


“k 


Jk 


* |uj 


|u:|At, 

2tt 


( 44 ) 


where |u£ is induced by outer flow minus the new born vortex. 
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Alternatively, this can be interpreted as a condition on the At^ 
of numerical integration of the boundary layer; i.e. 

2,nn k 


At. 


(45) 


uS 


'v 1/ ^ Re 2j)us * t ^ ie or ^ er t ^ ie boundary layer thickness 


where now 

(found, for example, by the first moment of the vorticity distribution). 

Now this argument must be compatible with the constraint on At^ 
imposed by the workability of the numerical integration scheme employed 
in the boundary layer. 

A At^ found by equation (45) was used successfully in early runs 

but it was found that increasing the At^ over this value lead to the 

same results with a savings in computer time. 

Thus the final results used a At. which was such that i) The 

k 

numerical integration scheme for the boundary layer equations was 
workable, and ii) the resulting m^, as given by equation (44), was 
"in" the boundary layer which was assumed to be much less than the 
radius of the cylinder; i.e. m^ < 0.1a. 


3.6 Rear Shear Layer 

The math ema tical model described so far leaves the satisfaction 
of the no-slip condition at the rear of the cylinder in the absence of 
the boundary layer (i.e., where the separation point is upstream) to 
the velocity induced by the array of free shed point vortices. These 
induce a velocity on the surface counter to that given by uniform 
flow over a cylinder. However, the counter-velocity sometimes 


exceeds the latter velocity and the no-slip condition is not 
satisfied. 

To satisfy the no-slip condition, the existence of a rear shear 
layer is postulated. The rear shear layer is similar to a boundary 
layer except that the outer flow is rotational. Since there is no 
theory for such a rear shear layer (time dependent) and since its 
prime contribution is the production of vorticity, then it is 
postulated: 

a) there exists a rear shear layer when |uo| > 0.1. 

For 0 < 6 < ir (with similar equations for tt < 0 < 2ir) , 

b) the separation point of rear shear layer is given as a 

function of the boundary layer separation point location 
by means of the following: (see Figure 6) 


0 - 6 


m 


e - e 
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e - e 

m s 


e - e 
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(46) 


Thus the location of the separation region of the rear shear 
layer is primarily a function of the boundary layer separation 
location. 

c) A point vortex is born at 9 of strength 
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u“ 2 At 

r r6k = F ’t "Vi— 
e=k 
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[a k + “rk’ 
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at a point 


(47) 


where 


l u o | At k 
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rk 2tt 


de 


only if 


sk 


dt 


d) N(< 5) point vortices are replaced by one point vortex of 
strength 


N 

r = I At 


k 2 


k=l 


at position m 


rk tt<uS> 


*■< 0.1 


< 0.1 


( 48 ) 


where 


<ul> = average of uj at 0 


This operation is used to reduce computer time and is 
allowable since one just born point vortex is of small 
strength (compared with one born from the boundary layer) 
and the separation point, as defined above, undergoes 
oscillation. Thus the operation is an averaging and co- 
alescing mechanism) • 

3 . 7 Wake-Vortex Convection 

The sum of point vortices in the outer flow define the wake. 
Each point vortex, labeled ft, moves with the cartesian velocity 
components of the outer flow. Thus 


x gk+l " X £k + At k u ^Bk’^k’V 
y 6 k+l = y 8k + At k ^Bk’^k’V 


( 49 ) 


where the contribution to [u°, v°] of the point vortex at 
is omitted. 
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3 « 8 Wake-Vortex Diffusion 


A point vortex undergoes both convection and diffusion. 
Convection is treated in the manner just described. Diffusion is 
handled as follows: for an isolated stationary point vortex, the 

governing equation is the diffusion equation with the potential 
solution as the initial condition. Thus following reference 39, 


DE 


2 r 

a). = 6 (to, + ) ; r > 0 

t .rr r 


t-0: r-0 


(50) 


0 , r > 0 


The solution is, with 6^ = l/Re(= l/Re 2DUg ) 


u,(r,t) = 4^t Re e 


-Re r /4t 


(51) 


with an induced velocity distribution (see Figure 7) . 


u(r,t) 


2 r 


2tt 


rr r _ 2 /» . 

U r dr d e " 2 ^ (1 - 6 } 


(52) 


The viscous core radius r (t) is defined as the distance from 

c 

the center to the point at which u is a maximum. From the above 
equations this is 


r (t) 

c 


V 


5.04t 

Re 


(53) 


For R - 4.10 


r (t) - 0.0112 Si 

c 


(54) 
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Then the potential flow solution for a point vortex is assumed 


valid for r > r c (see Figure 7) . 

The concept of a core radius r c is used in two applications. 

Surface pressure . If a point vortex is located in the neighborhood 
of the surface of the cylinder, a potential flow velocity distri- 
bution would induce a high velocity at the surface which would 
imply a peak in the pressure distribution on the surface. 

Approach of Two Point Vortices . If two vortices are close 
together anywhere in the flow, the potential flow velocity 
distribution would induce local high velocities. 

(Although peaks have been found experimentally, their small scale 
property is not compatible with the present mathematical model) . 

Both effects are resolved by the diffusion action which reduces 
the potential flow velocity in the neighborhood of the point vortex. 

One way to resolve this is to program an r^ as a function of the 
life-time of a point vortex via equation (53) and use the velocity 
distribution (see Figure 7). 

u = 0; r < r c (t) 

= —— ; r > r (t) (potential flow) (55) 
2trr c 

However, the employment of r^ as a function of t is insufficient 
within the present mathematical model which employs an approximation 
to the convection of the vortices as in equation (49) [i.e. Ax ^ At]. 
Thus the errors inherent in this convection approximation could place 
a point vortex sufficiently close to the surface or another point 



I 


vortex such that r for a small times is insufficient to rule out 

c 

high local velocities. 

In lieu of the r as a life-time function, the values of r * .05 
c c 

was used throughout. This is an average value of the life-time 
function, compatible with a large number of vortices, and is sufficient 
to rule out singularities. 


3.9 Wake-Vortex Coalescence 

The wake consists of a set of point vortices, the number of which 
grows with time. To reduce the computing time, the number is reduced 
at certain times by means of coalescing pairs of vortices into one 
vortex. The criteria are that the velocity field induced by the new 
vortex is equivalent to that induced by the original pair, and the 
times at which this operation is applied are those times at which the 
number of vortices is overly sufficient to yield an acceptable flow 
field and so a reduction of the number does not destroy the intrinsic 
structure of the wake (references 40 and 41) . 


The operation is as follows: two point vortices of strengths 

and r 2 , at [x^,y^] and l x 2> y 2^ respectively are replaced by one point 
vortex of strength at [x^y^] where 


(l r il*l +l r 2 l* 2 ) 
* 3 " (l r 1 l + l r 2 |) 

(l r il y i +lr 2 |y 2 > 

73 (1^1 + |r 2 j) 


(56) 
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(absolute signs to coalesce vortices of opposite sign). Then, if 

2 2 2 

(*2 ~ x i^ r i and ^ y 2 ” y i^ r i* r i ■ ( x “ x ± ) + (y - y ± ) = 1 » 2 * 

\ 

are of order e << 1 then 


u 3 » (u 1 + u 2 ) (1 + 0(e)) 

v 3 = (v 1 + v 2 ) (1 + 0(e)) (57) 

3*10 Wake-Vortex Stability 

For early times the flow field consists of an outer potential 
flow over a cylinder with free point vortices at the rear and an 
inner flow (the boundary layer and a possible rear shear layer) with 
all solutions symmetric about y = 0. It is known that this structure 
is unstable and asymmetry arises (references 42,43). This is usually 
obtained by supplying a small perturbation to the point vortex array 
and having the system determine the new stable asymmetric structure. 

However, for the times involved in the application of the theory, 
there would be negligible asymmetry. Thus to reduce computer time, 
symmetric wakes were assumed. 


3.11 Vortex Flux Parameter 

It was found that the accuracy of the resulting force and moment 
data in comparison with experimental data (see Section 5.0) could be 
improved in certain cases by reducing the vorticity in the wake by a 
factor (l-a)100% from the value actually produced by the separation 
of the boundary and rear shear layers. For the two tests performed 
(see Section 4.0) no reduction was required for the closed-ended 
ellipsoid (i.e. a = 1) but a reduction of 40% improved the results 
for the open-ended ogive -cylinder (i.e. a = 0.6). 
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This is the only empirical factor in the technique. Its role 
may be viewed as similar to the eddy viscosity in turbulence; i.e. 
empiricism arises in the absence of complete theoretical understanding. 
In this situation, there is no theoretical derivation of the 
equivalence underlying the unsteady cross flow analogy. Further 
light would be shed on this with further testing. 

This factor is not related to the ‘reduction in wake vorticity 
from that produced by the boundary layer alone in the classical 
problem of a circular cylinder of constant radius impulsively started 
from rest (see references 44,45). In earlier work with the constant- 
radius-solution, treating this problem, it was found that the rear 
shear layer effectively accounted for this reduction by introducing 
vorticity of sign opposite to that introduced by the boundary layer. 

3 . 12 Drag 

At representative times, the drag is computed in equations (11), 

(12) , (13) and (14) . The pressure is obtained from the Bernoulli 
equation, equation (35) and the shear stress from the boundary layer 
solution t = 26 u,- = 26 u(Ar,0t)/Ar by virtue of the no-slip condition. 

Within a completely inviscid model, the drag may be obtained using 
the vortex- impulse concept or Lagally's theorem. However, drag is 
obtained herein by pressure and shear integration since the wake 
model, while primarily inviscid, does have diffusive effects and, 
more importantly, the surface pressures are of interest in establishing 
the equivalence property. Computational times are of the same order. 



4.0 Applications 


4.1 Test Data 

The foregoing theory will be applied to two bodies for which 
experimental data is available. The first test case is that of an 
ogive-cylinder with a blunt base at subsonic speeds (reference 6) . 
The second test case is that of an ellipsoid of revolution at low 
speeds (reference 7). The latter presents a severe test in that the 
geometry is that of a closed-ended body such that for all angles of 
attack, there may not exist any boundary layer on the aft end. The 
theory as presented so far does not cover this eventuality (although 
it may easily be extended). However, even without this specific 
capability, meaningful results can be obtained. 


Test Case 1; Ogive-Cylinder (Reference 6) 

Geometry and Test Condition 

-. 1/2 
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; 0 < z* < a 
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- d (j^) -t]- d r(i- 
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z * 
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■ 4 » ^ t z* 1 1 

2 n 

where is the nose (ogive) length. 

£ - 1.282 m, d - .1194 m, l - 3d, f - 4 - 10.74 

n d 

Re 3DS - 4,7 x lo6, M " 0,3 
0 < a < 24° 

a - | - .0597 m 
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Invlscid Theory Results 


Slender Body 
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Test Case #2: Ellipsoid of Revolution (Reference 7) 

Geometry and Test Conditions 


r Q *(£*) = (|) v'zM* - z*) ; 0 < z* < £ 


% = 0.6096 m. , d - 0.1524 m. , f = ^ = 4 

d 

, . ~„2 , 
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v = 14.957 10 ^ m 2 /sec., Re^g = 4.97 x 10 5 
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Invlscid Theory Results 


Potential Flow 


c - f-f 0 cos 3 sin 2a 
n 1 l 


where where = Y ^2 ~ z ) > *2 = 
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16z(l - z) 
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Slender Body 


c - x (4 - z) sin 2a 
n l l 


C N “ °* C M,l/2 


2 sin 2a 
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4.2 Boundary Layer Integration Parameters 

The following parameter values were used: 

Ar = .14, A0 variable. At * .125 


4.3 Numerical Integration Parameters : 

The particular numerical integration sub-routines are to be 
found in Part II. The sub-routines are such that with a given integrand, 
less than 2% error is incurred in obtaining the integral. 


5.0 Results and Discussion 


5.1 General Results 

The results of applying the theory to the two test cases are 
to be found in Figures 8-18. 

Primarily the theory yields results capable of describing 
experimental data at high angles of attack where inviscid theories 
(potential and slender body) fail (Figures 8 and 14) . In achieving 
this state of data only one element of empiricism, the vortex flux 
parameter, a, was used. 

5 . 2 Ogive-Cylinder 

Figure 8 shows the basic results, the comparison of the values 
obtained by this technique and inviscid theory with experimental 



data for the normal force and pitching moment. The accuracy of 
inviscid theory, incapable of describing separation phenomena, 
decreases rapidly with increase in the angle of attack while the 
present theory remains valid. 

A value of 0.6 was used for the vortex flux parameter, a, to 
obtain these curves. There is the potential of increased accuracy 
in these results through further testing based upon variation of 
this parameter (see below, in discussion of Figure 10) . 

Figure 9 presents the basic input data to the curves of Figure 8, 

the sectional normal force coefficient c (z,cx) as a function of the 

n 

axial station z. It is noted that for 0 < z < 0.3, inviscid theory 
suffices. This behavior is pertinent to the discussion (see Section 2.1) 
of the property of no upstream influence intrinsic to the unsteady 
cross flow analogy. The flow field relevant to Figure 9 is subsonic 
which allows for upstream influence and the inviscid theory results 
contain this property. Yet there is agreement between the two, even 
at a = 20® where any slender body equivalence principle should not be 
expected to hold. This behavior lends support to the use of the 
unsteady cross flow analogy. 

For z > 0.3, inviscid theory predicts c ■ 0, while results of 

n 

employing the unsteady cross flow technique match experimental data. 

This is a separated region where the use of inviscid theory is 
obviously invalid. 

Figure 10 displays the effect of a change of a in the c (zyx) 
vs. z data. (At = .0625 was used in these runs while At ■ .125 
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was used in the runs of Figures 8 and 9. However, on the basis of 
early testing based upon At variation, the data is relatively 
insensitive to this change in At). Two values of a were used a = 0.5 
and a = 1.0. The reduction of a reduced the theoretical results 
to the neighborhood of the experimental data. An explanation for this 
behavior lies in the physics of the two dimensional unsteady flow 
problem. The drag value approaches zero as the wake disappears 
(D'Alambert's paradox). Thus a reduction in the wake strength 
reduces drag which is proportional to the sectional normal force. (Note 
that this argument, based on a cylinder of constant radius, is valid 
on the aft end of the ogive-cylinder). As a result of this test the 
value a = 0.6 was decided upon. 

It is noted that the experimental data was based upon a Re^g = 

Vi/v = 4.73 x 10 6 yielding the two-dimensional unsteady flow Reynolds 

4 

numbers R^pys = = 1*9, 3., 5.7, 7.5, and 8.9 x 10 for a = 5, 10, 

15, 20 and 24° respectively (see equation (16)). For a two dimensional 

/j, 

circular cylinder in steady flow a Reynolds number of 4.5 x 10 is used 
as a lower limit for transition and this value is exceeded for a 1 15°. 
However, the assumption is made that laminar boundary layers exist at 
all these Reynolds' numbers, including transition values. 

Figure 11 indicates typical behavior of the coefficient of 
pressure on the surface at various axial stations and angles of 
attack. As noted in the preceding the two-dimensional unsteady pressure 
distribution yields a logarithmic singularity. However, this term 
is independent of 0 and so a C , subject to an arbitrary additive 


1+2 



factor, can be determined. Thus these curves are qualitative. 

Within this context, the agreement is very good. 

Figure 12 presents the separation regions of three dimensional 
steady flow as obtained by this technique. In the absence of experimental 
data, these results can only be deemed physically acceptable. 

Figure 13 show wake vortex patterns. These present the positions 
of shed point vortices at one time or axial position and their position 
at the succeeding time or axial position. Thus the distribution of 
vorticity and the partial streamline pattern is indicated. 

5.3 Ellipsoid of Revolution 

The basic result for the ellipsoid, the comparison of the values 
obtained by this technique and inviscid theory with experimental data 
for the normal force and pitching moment about the midpoint are 
presented in Figure 14. In this case the improvement of this technique 
over inviscid theory is marked since the latter yields ■ 0 (reflecting 
the closed-end). 

These results employed a vortex flux factor of a ■ 1.0 (i.e. no 
reduction) (see discussion below for Figure 15) . 

The curves of Figure 14 are a result of the integration of the 
section normal force coefficient C n (z,o) over z. The integrand, 
c Q (z,a) is presented in Figure 15 as a function of the axial station £ 
at various angles of attack, a. In interpreting these curves as 
integrands, attention had to be paid to the behavior in the neighborhood 
of the base. In this region, there is no boundary layer; the tall 
lies within the wake. As noted earlier the theory as developed thus 
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far does not handle this situation although the extension can easily 
be made. And before the computations reach the section of no boundary 
layer, the radius a(t) may change rapidly (for a closed-ended body 
such as the ellipsoid) such that very fine integration meshes are 
required with consequent increase in computer time. Furthermore, the 
basic assumptions of the two-dimensional unsteady boundary layer 
solution technique may fail, particularly with regard to the property 
of no pressure change across the boundary layer. 

In lieu of increasing the complexity of the technique, computation 

was formally stopped at z = 0.9 □ z^ and at a straight line drawn to 

complete the curve at (c ,z) = (0,1.0), a theoretically exact point. 

n 

Figures 15b and 15c show the differences in behavior for two stopping 
stations z f = 0.9 and 1.0. The latter value, in both cases, departs 
from the experimental data. 

Figures 15b and 15c also show the effect of varying the vortex 
flux parameter a. For a = 10° there is little difference but for 
a = 20° the difference is marked. The behavior in this case is the 
same as for the ogive cylinder in that a reduction in a reduces drag 
which in turn reduces the section normal forces. However, in this 
case, the physics is not as clear since the corresponding two dimensions] 
unsteady problem is that concerned with a cylinder whose radius is 
changing rapidly with time. Also, the reduction in vortex flux 
Increases the error •’s opposed to the situation in the case of the 


3DS 


= 4.97 x 10' 


ogive-cylinder . 

The experimental data was obtained for a Re 



yielding the two dimensional unsteady flow Reynold's numbers of 
Re 2DUS = ^*0, ^.36 and 16.5 x 10^ at a - 4, 6, 10 and 20° 

respectively (see equation (16)). Based upon two dimensional steady 
flow, laminar boundary layers should exist. 

Figure l6 shows the pressure distribution for a * 20° at various z. 
The curve must still be considered as qualitative (i.e. with an 
arbitrary additive constant) . It is to be noted that qualitative 
agreement is found at z = .72 but for z > .86 where the tail is in 
the wake, even this qualitative agreement is lost. 

Figure 17 depicts the regions of separation of a three dimensional 
steady flow field as determined by this theory compared with that 
determined by a finite difference solution (reference 46) . The 
agreement is fairly good. Of note are the divergences on the aft 
windward region but in this region, the particular behavior determined 
by means of this technique agrees with experimental data (reference 47) . 

Figure 18 show wake vortex patterns, explained in the ogive- 
cylinder case. 

5.4 Vortex Flux Parameter 

As was noted, the value of the vortex flow parameter used was 
o = 1.0 for the ellipsoid and o ■ 0.6 for the ogive-cylinder. The 
prime difference between the two geometries is at the base: open- 
ended for the ogive-cylinder and closed-ended for the ellipsoid. 

Thus the recommendation for use of this technique for an arbitrary 
body is a value of a given by 


^5 



However, it is believed that substantial reductions are possible througl 
further work on the basic solution techniques and, to a lesser 
extent, in the programming. 

! 

6.0 Conclusions J 

The method developed herein employing the unsteady cross flow 
analogy and a discrete vorticity wake for the prediction of local flow 
conditions and overall forces and moments on bodies at high angles 
of attack with large regions of separated flow has been shown to be 
technically feasible. The method is based upon physical understanding 
of the flow and does not require a complect finite difference solution 
to the Navier-Stokes equations. 
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The advantages of such an approach lie in the potential of 
such solutions to bring about further physical understanding with 
consequent new analytic solutions (e.g. the wake) and in the area 
of computer-aided-design of aircraft where the engineer requires the 
physical understanding. In addition , there is the advantage of 
reduced computer time and storage although this may be offset somewhat 
with the advent of larger and faster computers. 

The results support the use of the unsteady cross flow analogy 
for three dimensional steady separated flows. A theoretically 
derived equivalence between two dimensional unsteady and three 
dimensional steady flow has been replaced by one empirical factor, 
a, with physical bounds, i.e. 0 < o £ 1 (cf« the eddy viscosity in 
turbulent flow). Further testing of the technique, with the use 
of equation (58) for a, should shed further light on the role played 


by this factor. 
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Figure 4 - Outer and Inner Flows 



Figure 5 - Comparison of Drag Coefficients for Impulsive 
Flow of cylinder of Constant Radius. 





Experiment (reference 6) 



Figure 8 - Aerodynamic Characteristics of Ogive-Cylinder 





Figure 8 - Concluded 




Experiment (reference 6) 




Figure 9 - Continued 
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Figure 9 - Concluded 







Present Theory (with arbitrary additive constant) 



Figure 11 - Continued 



Figure 11 - Continued 




Present Theory (with arbitrary additive constant) 




Figure 11 - Continued 


Present Theory (with arbitrary additive constant) 
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Figure 15 - Normal Force Distribution Over Ellipsoid 
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Figure 15 - Concluded 


Present Theory (with arbitrary additive constant) 



Figure 16 - Pressure Distribution on Ellipsoid 




Figure 16 - Continued 




Figure 16 - Concluded 
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Theoretical Separation Regions on Ellipsoid. 







